Gyroscopes based on nitrogen-vacancy centers in diamond 
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We propose solid-state gyroscopes based on ensembles of negatively charged nitrogen-vacancy 
(NV~) centers in diamond. In one scheme, rotation of the nitrogen- vacancy symmetry axis will 
induce Berry phase shifts in the NV~ electronic ground-state coherences proportional to the solid 
angle subtended by the symmetry axis. We estimate sensitivity in the range of 5 x 10"'^ rad/s/\/Hz 
in a 1 mm^ sensor volume using a simple Ramsey sequence. Incorporating dynamical decoupling to 
suppress dipolar relaxation may yield sensitivity at the level of 10~^ rad/s/v'Hz. With a modified 
Ramsey scheme. Berry phase shifts in the hyperfine sublevels would be employed. The projected 
sensitivity is in the range of 10~^ rad/s/\/Hz, however the smaller gyromagnetic ratio reduces 
sensitivity to magnetic-field noise by several orders of magnitude. Reaching 10~^ rad/s/\/Hz would 
represent an order of magnitude improvement over other compact, solid-state gyroscope technologies. 

PACS numbers: 03.65.Vf, 61.72.jn, 06.30.Gv 



Inertial sensors form the basis for navigation systems 
in both military and civilian applications. Such devices 
also find use in gravity research, for example, in tests of 
Einstein's equivalence principle [T] or to search for the 
Lens-Thirring effect [2], a prediction from general rela- 
tivity of "frame-dragging" in the presence of a massive 
rotating body. Present state-of-the-art sensors used for 
navigation in the airline industry are based on the Sagnac 
effect in fiber-optic bundles, with sensitivity at the level 
of 2 X 10~^ rad/s/-\/Hz j3j. The same effect in large-area 
(« m^) ring lasers has yielded sensitivity at the level of 
2 X 10^^° rad/s/\/Hz Gyroscopes based on cold atom 
interferometry have also demonstrated sensitivity in the 
range of 10"^ to 10"^° rad/s/VHz [213]. Noble gas nu- 
clear spins, which feature long coherence times, can also 
be used to form rotation sensors [61 [7], with demonstrated 
sensitivity of about 3 x 10~^ rad/s/-\/Hz. Variations on 
this scheme yield projected sensitivity in the range of 
-|^Q-io rad/s/v'Hz ^ in an active volume of several cm'^. 
This concept has been extended to miniaturized versions 
[S]. Finally, commercially available vibrating microelec- 
tromechanical systems (MEMS) gyroscopes achieve sen- 
sitivities on the order of 2 x 10~^ rad/s/-\/Hz [TUj. Ex- 
cept for the latter two cases, these sensors rely on large 
volumes or enclosed areas. 

Here we propose applying geometric-phase effects in 
diamond to form gyroscopic sensors from high density 
ensembles of NV~ centers. The key advantages of NV~ 
based sensors are their potentially small dimensions and 
thermal robustness, with the ability to operate from cryo- 
genic temperatures to 600-700 K [JJL . Recent theoretical 
work discusses the possibility of observing large geomet- 
ric phase shifts in single NV~ centers in diamond [T^ . 
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Similar ideas to those proposed here have been discussed 
in Ref. ^13^. Geometric phase shifts, also known as Berry 
phase shifts [T3] , arise when the symmetry axis of a sys- 
tem is rotated about an axis that is not parallel to this 
axis. This results in relative phase shifts of the mag- 
netic sublevels. One encounters Berry phase shifts, for 
example, when a spin system evolves in time in the pres- 
ence of an external magnetic field, which itself rotates in 
space in such a way that the field vector traces a nonzero 
solid angle [15]. In the case of adiabatic evolution, an 
eigenstate with magnetic quantum number m acquires a 
Berry phase (in addition to the phase due to evolution 
under the static Hamiltonian) that is given by the prod- 
uct of m and the solid angle traced by the field, dfl. If 
the quantization axis is determined by an electric field, 
a similar effect occurs on rotating the electric field [16) . 
This is the case for nuclear quadrupoles interacting with 
an electric-field gradient, and such an effect was used to 
demonstrate a ^'^^Xe nuclear-spin gyroscope [T7] . 

The NV^ center features a ground state with electronic 
spin S — 1 which can be initialized, manipulated, and de- 
tected via convenient optical and microwave transitions. 
Furthermore, the lifetime of the ground-state coherences 
can be quite long, on the order of 2 ms [TH] in ^'^C de- 
pleted samples, even without applying dynamic decou- 
pling sequences. The analysis presented below indicates 
that sensitivities at the level of ~ lO"'^ rad/s/-\/Hz should 
be achievable in mm scale packages. This falls short of 
the sensitivities demonstrated by state-of-the-art tech- 
niques based on larger sensors, though it surpasses that 
of commercially available vibrating MEMS sensors. Fur- 
thermore, operation as a three-axis sensor may be pos- 
sible by monitoring coherences in NV^ centers with dif- 
ferent orientations with respect to the host lattice. 

The ground-state level structure of the NV~ centers is 
shown in Fig. [ija). Detailed discussions of the physics 
and spectroscopy of this defect can be found in Refs. 
[TUH^ and references therein. NV~ centers possess C^y 
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symmetry, and can be oriented along any of the four crys- 
tallographic orientations. The electronic ground state is 
spin-triplet (S = I). In the absence of strain or exter- 
nal magnetic field, the ms = ±1 levels are degenerate, 
where ms is the quantum number associated with the 
projection of S along the symmetry axis. Due to the 
spin-spin interaction, the zero-field splitting between the 
ms — and ms = ±1 levels is D = 2.87 GHz. Appli- 
cation of a magnetic field along the NV~ symmetry 
axis produces a shift of the ground-state energy levels ac- 
cording to msgsfJ-sBz, where « 2 is the Lande factor 
for the electron and — 8.79 rad/s/G is the Bohr mag- 
neton. (We focus here on a single orientation, assumed 
parallel to the magnetic field. The projection of such a 
magnetic field onto the other three axes is suppressed by 
a factor of 3. Therefore, for sufficiently large magnetic 
fields, ~ 10 G, the microwave transitions of differently 
oriented NV^ centers are well isolated from NV~ cen- 
ters whose symmetry axis is parallel to the field. Since 
microwave transitions are excited resonantly, we can ig- 
nore the contributions from these other orientations.) In 
the presence of local strain, the mg = ±1 sublevels are 
mixed and the corresponding eigenstates are split in en- 
ergy by 2E, where E is the transverse zero-field splitting 
parameter, typically on the order of a few MHz for high 
NV~ density samples. Application of magnetic fields of 
several gauss will overwhelm this, so that eigenstates are 
approximately those of Sz ■ Finally, hyperfine coupling of 
the ^^N nuclear spin (/ — 1) to the electronic degrees of 
freedom produces additional splittings as shown at the 
rig ht of Fig. [TJa). 

Optically detected magnetic resonance (ODMR) spec- 
troscopy of the ground state can be accomplished as fol- 
lows: Optical transitions from the ground states to 
the excited-state manifold are spin-conserving, pre- 
serving the magnetic quantum number ms, and are typ- 
ically excited via the phonon sideband (PSB) of the 
'^A2 — J-^ E transition using 532 nm light. Optical pump- 
ing of the ms = ground state occurs through spin 
non-conserving decay pathways (due to spin-orbit cou- 
pling) involving singlet levels not shown in Fig. [T] Flu- 
orescence in the zero-phonon line (ZPL) at 637 nm and 
the associated PSB is brighter for NV~ starting in the 
ms = state. Thus, application of microwaves tuned 
to the ms = 0— >m5 = ±l ground-state transitions 
depopulates the ms = state and reduces fluorescence 
in the ZPL and its PSB. Assuming the NV~ center is 
aligned with the magnetic field, and neglecting ^"^N hy- 
perfine coupling, fluorescence as a function of microwave 
frequency appears as a doublet, centered about D, with a 
splitting determined by the magnetic field A — 2gsiiBBz 
(assuming QsIJ-bBz ^ E). If ^''N hyperfine coupling is 
resolved, each line of the doublet splits into a triplet. 
ODMR spectra are schematically shown in Fig. [ijb) for 
unpolarized and polarized ^^N spins (top and bottom 
traces, respectively). The polarization of nuclear spins 
involves elevated magnetic fields (see below), however, 
this plot assumes observation in a magnetic field w 10 G 
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FIG. 1: (a) Ground-state energy levels of NV~ centers in the 
presence of a magnetic field of several G, and accounting for 
^^N hyperfine coupling. The solid, double-headed red arrows 
show the allowed microwave transitions, which preserve the 
nitrogen spin. Ami = 0. (b) ^''N hyperfine resolved ODMR 
spectrum as a function of microwave detuning from D for the 
case of unpolarized (upper trace) and polarized (lower trace) 
^''N spins. 



in both cases. One can also detect the spin state by mon- 
itoring transmission of infrared light tuned to the 1042 
nm transition between the two singlet levels, resulting in 
improved light collection [22 . 

NV~ electronic Berry phase based gyroscopes - We first 
neglect "'^^N hyperfine coupling and suppose an ensemble 
of NV~ centers is pumped into the ms — state, with a 
magnetic field (approx 10 G) applied along the NV~ axis. 
Suppose the diamond is rotated with an angular velocity 
Lo about an axis z' oriented at an angle relative to the 
NV^ symmetry axis {0 < 7r/2) during some interval t 
(Fig. [2]). The symmetry axis thus traces a solid angle 
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d<j>d0'sm0' =ujt{l-cos0). (1) 
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In the course of this rotation, the Zeeman eigenstates of 
the NV~ center pick up a geometric phase shift A(f> — 
m^ri. From Eq. Q, we see that this geometric phase 
shift is maximized for 6* = 7r/2. If one establishes co- 
herences between ms — and ms — ±1 states prior 
to rotation, the observed transition frequencies will be 
u!± — 2ttD ± gsfiBB ± w(l — cos 0). 

The observed frequency shifts associated with rota- 
tions involves both the rotation rate uj and the angle of 
rotation with respect to the NV~ symmetry axis, 0. If 
the magnetic field is oriented such that resonances due to 
all four NV~ orientations are resolved. Berry phase shifts 
in the four different orientations can be used to determine 
both the direction of rotation and its magnitude. 

One can envision several schemes for readout. Here, 
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FIG. 2: Rotation of diamond about an axis z' oriented at an 
angle 9 with respect to the NV~ symmetry axis. 



we consider a Ramsey scheme, in which NV^ centers are 
initiahzed in the ms — sublevel by apphcation of 532 
nm pump Ught. A tt/2 microwave pulse is applied to the 
spins, they are allowed to precess freely for a time t, and 
then a second 7r/2 pulse is applied, converting coherences 
back into populations. Finally the spin state is read out 
by monitoring fluorescence while illuminating the sample 
with 532 nm pump light. The resulting fluorescence is 
given by 



F rjN {1 + R cos uj±t). 



(2) 



Here, R is the detection contrast for states with ms = 
and ms = ±1, rj is the collection efficiency, and N is 
the number of participating NV^ centers. This assumes 
a roughly optimal readout, in which each center emits 
one photon in the readout stage. The slope of the flu- 
orescence with respect to rotation rate lu (assuming lu 
is small) is maximized by tuning B and t such that 
{2ttD + gsfJ.BB)t — {n zL l/4)27r, where n is an inte- 
ger, thereby optimizing the sensitivity to rotations. Un- 
der such conditions dF/duj « ztrjNRt. The accuracy 
with which one can determine the rotation frequency w 
is 5uj = {dF/dLd)~^SF. Fundamentally limiting the noise 
in the fluorescence is photon shot noise, SF — y/rjN. Set- 
ting t — T2 , and repeating the measurement r/Tj times 
(t is the total integration time), we find 
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Realistic parameters may be conservatively estimated 
based on Ref. 23 , where NV~ ensembles were realized 
with density oi n — IQ^^ cm^'^, « 300 ns, and a 
measurement volume of 1 mm'^. Fluorescence detection 
suffers from poor contrast ratio and poor detection effi- 
ciency (typically R w 0.03 and 77 « 0.01). The detec- 
tion efficiency can be substantially improved to 1 
by employing IR absorption [22] . or by using the side 



collection technique demonstrated in Ref. [21], where 
r] ss 0.5 was realized. Assuming 77 = 0.5, 1/4 of the NV~ 
centers participating, a measurement volume of 1 mm'^, 
and a measurement time r = 0.5 s (corresponding to 
a bandwidth of 1 Hz), Eq. ^ yields rotation sensitiv- 
ity of about 5.4 x 10^'^ rad/s/\/Hz. This is below the 
demonstrated sensitivity of the state-of-the-art sensors 
mentioned above, however, the active part of the sensor 
is smaller by several orders of magnitude. 

The numbers given above are fairly conservative. Ref- 
erence [25] reports that CPMG decouphng sequences can 
be used to extend coherence times by an order of magni- 
tude when ^^N impurities limit the lifetime. More recent 
work has shown coherence lifetimes as long as 300 ms 
in low NV~ density samples in isotopically pure ^^C di- 
amond using CPMG decoupling sequence at 77 K [53]. 
For the work reported in Ref. [35] , dipole-dipole coupling 
between NV~ centers appeared to limit the lifetime. As 
suggested in Ref. [27], such relaxation may be mitigated 
by pulse sequences such as MREV [55] or WAHUHA [25] . 
which average the dipole-dipole interaction to zero and 
scale the rank-1 gyroscopic phase shift (for WAHUHA, 
the scaling factor is 1/ a/3) . Ultimately, Ti places an up- 
per bound on T2. In Ref. ^30J room-temperature mea- 
surements of Ti = 3 ms was measured for the high NV~ 
density samples of present interest. Assuming dipole- 
dipole relaxation can be sufficiently suppressed in high- 
density ensembles, one can expect roughly two orders 
of magnitude improvement in sensitivity over the esti- 
mate given above, placing the sensitivity in the range of 
10-5 rad/s/A/Hz. 

Magnetic-field fluctuations present a serious source of 
noise in this scheme. The sensitivity to magnetic field 
fluctuations is given by = {dF/duj)^^{dF/dB)SB = 
gsfJ-sSB. In order to reach an angular- frequency sensi- 
tivity of 5.4 X 10"'^ rad s"^, one requires magnetic field 
stability SB ss 3 x 10^^° G. Assuming w 10 G is required 
to sufficiently suppress the stress parameter, this corre- 
sponds to fractional field stability at the level of 10"^°, 
challenging, though not insurmountable. 

The temperature dependence of the D coefficient 
places substantial requirements on thermal stability. The 
temperature dependence of D near room temperature is 
linear [31], with a slope of about -70 kHz/K. To reach 
10^^ rad/s/-\/Hz requires thermal stability at the level of 
10"^ K. Such stability seems unrealistic, however, if one 
monitors the splitting of the ODMR signal for ms = ±1 
rather than the position of just a single line, the tem- 
perature dependence should be eliminated. This also im- 
proves the sensitivity of the gyroscope by a factor of \/2. 

Nuclear spin based diamond gyroscopes - Nuclear spins 
are appealing in the context of gyroscopes because their 
gyromagnetic ratio is about three orders of magnitude 
smaller than that of the electron, relaxing the require- 
ments on field stability. Recent work has shown that it is 
possible to achieve very high, « 98%, polarization of sin- 
gle ^^N or ^^N nuclear spins in the NV" center (as well as 
nearby ^'^C spins) using an excited-state level anticross- 
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FIG. 3: Pulse sequence for Berry phase based gyroscope using 
resolved '^''N hyperfine structure in diamond. 



ing at about 500 G, induced by the hyperfine coupling of 
the nucleus with the electrons of the center . This 

technique has been extended to polarization of ^^N and 
proximal ^'^C nuclear spins in ensembles of NV~ centers 
[33] . The ^'*N spin state can be detected via fluorescence 
or IR absorption by the NV^ center. 

An experimental protocol for a ^^N diamond based gy- 
roscope is presented in Fig. [s] NV~ centers are polarized 
in the {ms = 0,m/ — —1) state by turning the magnetic 
field (oriented parallel to the NV~ symmetry axis) up to 
about 500 G, while illuminating the sample with 532 nm 
light (see Refs. [SUM]). The hght is turned off, and the 
magnetic field reduced to about 10 G (to improve sta- 
bility), and coherences between the \ms = 0, m/ — —1) 
and \ms — 0,mi — 0) ground state sublevels are estab- 
lished by application of a tt/2 pulse, resonantly tuned 
to the hyperfine transition frequency ujhf ~ 5.1 MHz. 
Spins are allowed to evolve freely for some time t in the 
presence of the magnetic field, hyperfine coupling, and 
sample rotation, accruing both nuclear Zeeman, hyper- 
fine, and geometric phase shifts. A second 7r/2 pulse 
converts these coherences back into ^''N Zeeman-sublevel 
populations, which are then selectively read out by ap- 
plying microwaves tuned to the desired ^^N hyperfine 
transition, indicated by either of the downwards point- 
ing arrows in Fig. [ijb). Population in the selected level 
would oscillate sinusoidally with a frequency proportional 
to LOnuc = ^hf + "/nB + Lo whcrc is the gyromagnctic 
ratio of the nitrogen nucleus (we have set 9 = 7r/2 so that 
the factor 1 — cos 9 in Eq. ([T]) is equal to 1). Since the mi- 
crowave and optical transitions preserve rn/, this results 
in a sinusoidal modulation of the resulting fluorescence 
F 



Here, Bq is the magnetic field corresponding to resonance 
with the applied microwaves and Ai? = l/(r25sMs) is 
the width of the microwave transition in magnetic field 
units. The Lorentzian profile for the microwaves is in- 
cluded to address the possibility of noise due to magnetic 
field fluctuations. As we see below, it turns out that this 
factor drops out. The same arguments used to derive the 
sensitivity, Eq. ([s]), of gyroscopes based on pure elec- 
tronic transitions applies in this case as well. 



The exact value of will depend on conditions and 
sample preparation. We point out that in some instances, 
the nuclear spin coherence times r2_„„c can be signifi- 
cantly longer for nuclear spins than for electron spins. 
For example, with single NV~ centers, Ref. [35] reports, 
Tl^ « 1.4 /^s, while r2*n«c « 58 for spins. In 
^^C depleted diamonds with low NV~ density, ^^N coher- 
ence lifetimes T2 = 7 ms have been realized [36]. Mea- 
surements of these quantities in high-density samples are 
planned. Taking TJ„„c — 1 '^s, equation ([s]) yields sen- 
sitivity at the level of 9 x 10^^ rad/s/\/Hz, with other 
parameters as used above, and without requiring dynam- 
ical decoupling. 



F (xl + R 
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(4) 



Despite the fact that microwave transitions are used 
to selectively read out the nuclear spin state, the sen- 
sitivity to magnetic field fluctuations is dramatically re- 
duced compared to the case of gyroscopes based solely on 
electronic transitions: 5uj = {dF/duj)~^{dF/dB)SB = 
'YnSB. We see that this is suppressed relative to the case 
of electronic coherences by a factor of Jn /{gsf^B) ~ 9000 
for ^*N, greatly relaxing the requirements on field stabil- 
ity. 



In summary, we have proposed a new type of solid- 
state gyroscope based on Berry-phase effects in diamond 
nitrogen-vacancy centers. We estimate that the sensi- 
tivity of such a device could surpass that of commercial 
vibrating MEMS sensors, which find wide application in 
a variety of technologies. We also point out that NV~ 
centers have been used to polarize bulk ^^C in elevated 
magnetic fields [33 , and efforts are under way to pro- 
duce polarization of bulk ^^C nuclear spins in low mag- 
netic field. Bulk ^"^C spins can potentially have very long 
Ti times [35], and if decoupling schemes can yield cor- 
respondingly long T2 times, one may expect dramatic 
improvements in sensitivity. 
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